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Abstract. Submillimeter observations are a key for answering many of the big questions in 
modern-day astrophysics, such as how stars and planets form, how galaxies evolve, and how ma- 
O . terial cycles through stars and the interstellar medium. With the upcoming large submillimeter 
facilities ALMA and Herschel a new window will open to study these questions. ARTIST is a 
C/^ ' project funded in context of the European ASTRONET program with the aim of developing a 
^ I , next generation model suite for comprehensive multi-dimensional radiative transfer calculations 
of the dust and line emission, as well as their polarization, to help interpret observations with 
these groundbreaking facilities. 
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1. Introduction 

The Atacama Large Millimeter Array is the largest ground based project in astron- 
^-H ' omy. It will be the world's most powerful instrument for millimeter and submillimeter 
T-H I astronomy, providing enormous improvements in sensitivity, resolution and imaging fi- 
^ • delity in these wavelength bands. The main focus of the early use of ALMA, compared 
•'^ ' to present-day facilities, will be its high angular resolution and sensitivity. For studies 
/\ of star formation, for example, ALMA will zoom in to AU scales in circumstellar disks 
^ in nearby star forming regions and thereby address some of the key questions in disk 
" " ■ formation. It will also open up new possibilities in the study of planet formation, re- 
solving the effects of planets on the disks around young stars and thus providing direct 
observational constraints on planet formation models. At the same time, the Herschel 
Space Observatory is providing high spatial and spectral resolution observations at wave- 
lengths unobservable from the ground, in particular, of lines of II2O and high excitation 
transitions of molecules tracing the dense warm gas, e.g., in the innermost regions of 
protostars or in shocks. 
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The simultaneous observation of dust and of a multitude of spectral lines opens new 
horizons for the physical and chemical analysis of the objects, e.g., to determine the 
excitation conditions (density, temperature, radiation), the chemical abundances and 
chemical network, and through polarization measurements, the magnetic field. It will be 
necessary to take all these observational constraints into account for a realistic quantita- 
tive description. ALMA will offer a new chance to study magnetic fields: due to receiver 
constraints, full polarization calibration and imaging will be the norm rather than the 
exception, making it essential that a self-consistent polarization modeling tool is available 
to the ALMA users. 

With the novelty of these observations it will be critical to have an efficient, flexible 
and state-of-the-art modeling package that can provide a direct link between the theo- 
retical predictions and the quantitative constraints from the submillimeter observations. 
The new observational opportunities require a new generation of modeling tools that 
can model the full multi-dimensional structure of, e.g., a low-mass protostar, including 
its envelope, disk, outflow and magnetic field, and their time evolution. Current tools 
are inadequate for the modeling of such complex structures because of their speed and 
inaccessibility, while tools to model polarization are completely lacking. Both ALMA and 
Herschel will provide us with large samples of sources, observed homogeneously as part 
of large key and legacy programs. This makes it prudent to have easily accessible and effi- 
cient tools, which with high convergence speed incorporate all observational constraints, 
for large source samples, and in a systematic fashion. It is the aim of this program to 
provide such an innovative suite of model tools and test it with existing submillimeter 
data and with new data from ALMA and Herschel. 

2. Objectives 

The goal of this project is to deliver a next generation radiative transfer modeling 
package that provides a self-consistent model for the emission of a multi-dimensional 
source observed at submillimeter wavelengths. For this project we are specifically moti- 
vated, without any loss of more general applicability, by low-mass star formation. Our 
modeling package shall be able to provide a self-consistent modeling tool for the line and 
continuum as well as polarization emission from, e.g., a young stellar object, incorporat- 
ing an infalling large-scale envelope, rotationally decoupled protoplanetary disk, outflow 
cavity, and magnetic field - with no restrictions of the intrinsic geometry. With such 
an innovative model tool it will be possible to provide quantitative constraints on the 
relation between the large-scale angular momentum in the core and the disk evolution, 
on the direct impact of the outflows and their launching close to the disk surface, and on 
the importance of the magnetic field. Two important issues need to be addressed in this 
particular example: (i) young stellar objects are characterized by structure on a wide 
range of spatial scales and with complex geometries, but current radiative transfer tools 
are locked to fixed linear or logarithmic grids and can therefore model multi-dimensional 
source structures only with great computational time expense; (ii) while most current 
observations indicate that magnetic fields play an important role in various stages of 
the star formation process, their relation to the physical source structure is yet poorly 
constrained. ALMA observations will study the magnetic fields through resolved line 
and continuum polarization observations that can not be properly analyzed with current 
models. The modeling tool will also be useful for tackling scientific questions relating to 
ALMA observations of, e.g., evolved stars, planetary nebulae or extragalactic starbursts. 
Providing radiative transfer tools that take complex source structures into account will 
also be of great help to interpret Herschel observations, e.g., to understand the origin of 
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H2O and high excitation CO hnes in the interface between protostehar cores and outflows 
or jets. There is a high degree of couphng between different modehng aspects: e.g., to 
understand the polarization of a given molecule's emission it is necessary to understand 
the physical conditions and chemical network that leads to the molecule's formation 
and excitation. For complex source structures it is necessary to develop an approach to 
time-dependent multi-dimensional modeling, e.g., through libraries of theoretical model 
prescriptions that can readily be incorporated and expanded for comparison to the data. 



3. Tools 

The planned model suite will have the following three components: (a) an innovative 
radiative transfer code using adaptive gridding that allows simulations of sources with 
arbitrary multi-dimensional and time-dependent structures ensuring a rapid convergence 
and thus allowing an exploration of parameters; (&) unique tools for modeling the polar- 
ization of the line and dust emission, information that will come with standard ALMA 
observations; (c) a comprehensive Python-based interface connecting these packages, thus 
with direct link to, e.g., ALMA data reduction software (CAS A). A schematic overview 
of the program is shown in Fig. [TJ 

(a) — 3D Line Radiative Transfer: LIME. ALMA's high resolution data will produce 
the need to model phenomena with non-symmetric structures, such as spiral-waves, 
proto-planet resonances in evolving circumstellar disks, close protostehar binaries etc. 
Conventional radiative transfer tools use simple linear or logarithmic spatial grids. To 
model complex source structures in higher dimensions thus requires increasingly finer 
grid scaling, which becomes very difficult to handle computationally. As an alternative, 
the Simplex algorithm (Ritzerveld & Icke 2006) uses a Poisson method to define a grid 
based on the density distribution. The cornerstone of ARTIST is the 3D line radiative 
transfer code, LIME (Brinch & Hogerheijde 2010), which utilizes the SimpleX gridding 
algorithm in a 3D extension to the RATRAN radiative transfer code (Hogerheijde & van 
der Tak 2000). LIME is currently being applied to for example modeling new H2O obser- 
vations of protostars and disks from the Herschel Space Observatory. For these models 
the adaptive gridding method ensures rapid convergence, for lines from molecules such 
as II2O that are far from LTE excitation. 

(h) — Dust and Line Polarization. Various theoretical studies predict that magnetic 
fields, turbulence or/and magneto-hydrodynamic waves may be the main agents control- 
ling both the evolution of molecular clouds and the star-formation process (e.g., Bertoldi 
& McKee 1992; Mac Low & Klessen 2004; Mouschovias et al. 2006; van Loo et al. 2007). 
Unfortunately, the magnetic field is the least-known observable in star formation, due to 
the inherent difficulty to measure it with present telescopes (mainly through polarimetric 
observations of dust and molecular emission; e.g., Girart, Rao & Marrone 2006). This 
situation will change dramatically with ALMA, which will provide such an improvement 
in sensitivity that polarization observations of dust and molecular lines can be done in 
many sources at a very good angular resolution. However, the interpretation of the polar- 
ized data is difficult and appropriate tools are needed to scientifically harvest the wealth 
of polarization data expected from ALMA. Note that the ALMA band-7 (275 - 373 GHz) 
receivers will require full polarization calibration and imaging even for projects with no 
primary interest in polarization. The ARTIST package consists of a set of modeling tools 
for polarization and magnetic fields in different molecular and dust environments (e.g., 
low- and high-mass protostars, envelopes around evolved stars). 
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Theoretical input models 

e.g., analytical collapse, simulation output (incl. magnetic field, chemical network) 



Dust radiative transfer fextema/ioo/sj 

Self-consistent (dust) temperature profile, dust density distribution 



Line excitation (3D) 

Chemistry: abundances 
Dynamics: velocity field 



Continuum polarization 

Grain alignment efficiency 
dependent on dens. / temp. 



Line polarization 

magnetic field model 
Zeeman splitting 



Ray-tracer 

images of molecular lines, continuum and polarization 




Integration of infiages into data analysis software, e.g., CASA 

Constraints on: density, temperature, velocity field, magnetic field, chemistry 



Figure 1. Schematic overview of the components in the ARTIST program. 

(c) — Model interface and library. An important component of ARTIST is a common 
interface for the codes used for radiative transfer modeling of typical data produced by 
submillimeter telescopes such as ALMA and Herschel. The model interface will include: 
(i) a library of standard input models, e.g., for collapsing protostars, circumstellar disks 
and evolved stars; (m) a wrapper package that links the existing dust and line radiative 
transfer codes; (Hi) tools to analyze model output, e.g., to extract information about 
molecular excitation and its deviation from LTE, or optical depth surfaces, and to import 
this in existing visualization packages; (iv) a ray tracing backend that readily provides 
data cubes that can be used in data reduction packages. 



4. Current Status 

The aim of the ARTIST project is to supply the community with the described tools 
in one coherent modeling package. The tools will be made publicly available as they are 
finished: we expect LIME to be released medio-2010 with the remaining components of 
the package to follow. The tools will be distributed and supported through the ALMA 
regional center nodes in Bonn and Leiden as well as the Danish initiative for Far-infrared 
and Submillimeter Astronomy (DFSA; Copenhagen, Denmark). For more information 
see http : // www . astro . uni-bonn . de/ ARC/ artist . 
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